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Abstract 

Using the data coming from the new 182 Gold type la supernova samples, the shift parameter of the Cosmic 
Microwave Background given by the three- year Wilkinson Microwave Anisotropy Probe observations, and 
the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, H{z) and lookback time 
measurements, we have performed a statistical joint analysis of the interacting holographic dark energy 
model. Consistent parameter estimations show us that the interacting holographic dark energy model is a 
viable candidate to explain the observed acceleration of our universe. 
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The present continuous flow of cosmological data provides us day by day a clearer picture that 
our universe is experiencing accelerated expansion [1]. In order to draw precise conclusions from 
the available phenomenological perspectives on how fast the universe is expanding at present, how 
long this speed up has lasted and how the acceleration rate has changed over the recent past, we 
have to overcome statistical uncertainties and possible theoretical biasing in the tests used. This 
requires us to refine the existing tests and devise new ones. 

Recently Simon et al have published Hubble parameter data extracted from differential ages 
of passively evolving galaxies. It is interesting to use these data to constrain the evolution of 
the universe. This is so because that they can provide consistent checks and tight constraints on 
models when combined with other cosmological tests, and also because the Hubble parameter is 
not integrated over like that of the luminosity distance and it can give better constraints on the 
cosmological parameters. Recently, Hubble parameter data have been used to constrain several 
cosmological models [3j, |4 1 . 

To reduce the degeneracy in viable candidate cosmological models designed to explain the 
observed accelerated expansion, new observables should be added to the usual ones. Recalling that 
the test of cosmological models by the type la supernova (SN la) data is a distance based method, 
it is of interest to look for tests based on time-dependent observable. In 0,0], the age of an old high 
redshift galaxy has been used to constrain the model. To overcome the problem that the estimate 
of the age of a single galaxy maybe affected by systematic errors, it is needed to consider a sample 
of galaxies belonging to the same cluster. Recently, the age estimates of around 160 galaxy clusters 
at six redshifts distributed in the interval 0.10 < z < 1.27 have been compiled by Capozziello et 
al Q]. Employing these data, one can take into account the lookback time which was defined by 
Sandage [s] as the difference between the present age of the universe and its age when a particular 
light ray at redshift z was emitted. This quantity can discriminate among different cosmological 
models. The lookback time has been used as a test for some cosmological models [?, ^. 

In this paper we will use ..e .a.e. SN .a d.a compiled Riess e. a> Q. .he CoL. Micowave 
Background (CMB) shift parameter derived from the three-year Wilkinson Microwave Anisotropy 
Probe (WMAP3) observations Ul], the baryon acoustic oscillations (BAO) measurement from the 
large-scale correlation function of the Sloan Digital Sky Survey (SDSS) luminous red galaxies 12 1 
in combination with the H(z) data and the lookback time data to give a coniplete investigation 
on the viability of the interacting holographic dark energy model devised in [13^. Recently, this 
model has confronted the tests from the SN la data [l^, the age constraint and the small / CMB 
spectrum constraint . It has been argued that the interacting holographic dark energy model can 
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accommodate the transition of the dark energy equation of state w from w > —Itow < — 1 [l^, [3] j 



as recently revealed from extensive data analysis [15|. With the interaction between dark energy 
and dark matter introduced in [l^, it has been shown that the old astrophysical structures can be 
formed naturally [6] and the coincidence problem can be alleviated 6j, [la] • The thermodynamical 
properties of the universe with the interacting holographic dark energy have also been studied • 
Very recently, the combined constraint on the interacting holographic dark energy model using 
the SN la data, the BAO measurement and the shift parameter determined from the SDSS and 
WMAP3 has been reported [3]. This paper aims to place combined new observational constraints 
on this interacting holographic dark energy model by including the Hubble parameter data and 
the lookback time data. Different from the distance based test, the lookback time is a time based 
method. Moreover, the Hubble parameter does not suffer the integration effect in the luminosity 
distance. It is expected that these new tests will further constrain the model. 

Recently, inspired by the holographic hypothesis [l^, a new model has been put forward to 
explain the dark energy. The energy density cannot exceed the mass of a black hole with the 
same size of the universe L, thus we have pD = 3c^L~^, where c is a constant and the Planck 
mass Mp has been taken unity. Choosing L as the future event horizon, Rfi = a "^J) have 
Pd = 3c^Rj^^ as the dark energy density. As far as energy conservation is concerned, we suppose 
that the interaction is described by the (separately non conserving) equations 

Pm + 3Hpm = Q (1) 

ph + SH{l + u;D) = -Q (2) 

where Q is some interaction term. For the moment we take for granted that the interaction is the 
one proposed on general grounds in [20], which is Q = 3b'^H{pm + Pd), where 6^ is the second 
phenomenological constant indicating coupling between dark energy and dark matter. Positive 
values of 6^ would correspond to a transfer of energy from the dark energy to dark matter, while 
the negative would imply a transfer of energy from the dark matter to the dark energy [211]. In 
view of the unknown nature of dark matter and dark energy, we do not put any limit on the sign of 
6^ at first and wait to determine it from the observational data. Because of the interaction, neither 
dark energy nor dark matter conserve whence they evolve separately. For the flat universe, using 
the Priedmann equation + Q^n = 1) where = and Q.m = ^2-, the evolution behavior of 
the dark energy was obtained as jl^ : 
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Figure 1: (a)The contours from the combination of SN la, BAO in the interacting holographic dark energy 
model for c and Q,dq at Icr, 2<j, confidence level with — —0.10. (6)The contours from the combination 
of SN la, BAO for and VIdq at Icr, 2(T, Sct confidence level with c = 0.53. 



The prime denotes the derivative with respect to x = In a. The equation of state of dark energy 
was expressed as 

1 



3c 



(4) 



By suitably choosing the coupling between dark energy and dark matter, this model can ac- 
commodate the transition of the dark energy equation of state from ujo > — ItowD < —1 jisl, [l^. 
which is in agreement with the recent analysis of the SN la data [15]. The deceleration parameter 
has the form 
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The evolution of the Hubble parameter can be written as 

l + q' 



Hz) 



Hq exp 



dz'] 



(6) 



/o 1 + ^' 

Next we constrain the interacting holographic dark energy model by using the latest observa- 
tional data, such as the gold SN la data, the shift parameter and the BAO measurement from 
WMAP3 and SDSS, and combining these observations with H{z) data and lookback time data. 
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Figure 2: (a)The contours from the combination of SN la, BAO, CMB in the interacting holographic dark 
energy model for c and JIdq at la, 2a, 3a confidence level with = —0.004. (&)The contours from the 
combination of SN la, BAO, CMB for and floo at la, 2a, 3cr confidence level with c — 0.84. 
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Figure 3: (a)The contours from the combination of SN la, BAO, H{z) in the interacting holographic dark 
energy model for c and JIdq at la, 2a, 3a confidence level with b^ = —0.005. (6)The contours from the 
combination of SN la, BAO, H{z) for b^ and JIdq at la, 2a, 3a confidence level with c — 0.82. We have 
employed Hq = 12km ■ s^^ ■ Mpc^^ . 
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Figure 4: (a)The contours from the combination of SN la, BAO, H{z), CMB in the interacting holographic 
dark energy model for c and fi/^o at la, 2a, 3a confidence level with = —0.003. (6)The contours from 
the combination of SN la ,BAO, H{z), CMB for and Q.do at Ic, 2(T, 'ia confidence level with c = 0.84. 
We have employed ifo = 72km ■ • Mpc~^. 
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Figure 5: (a)The contours from the combination of SN la, BAO, Lookback time in the interacting holographic 
dark energy model for c and JIdo at la, 2a, 3a confidence level with = —0.059. (6)The contours from the 
combination of SN la, BAO, Lookback time for and flnn at la, 2a, 3a confidence level with c = 0.62. 
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Figure 6: (a)The contours from the combination of SN la, BAO, Lookback time, CMB in the interacting 
holographic dark energy model for c and at Icr, 2a, 3(t confidence level with fe^ ~ —0.003. (6)The 
contours from the combination of SN la, BAO, Lookback time, CMB for 6^ and fi^jQ at la, 2a, 3a confidence 
level with c — 0.83. 

The up-to-date gold SN la sample was compiled by Riess et al [l^. This sample consists of 
182 data, in which 16 points with 0.46 < z < 1.39 were obtained recently by the Hubble Space 
Telescope (HST), 47 points with 0.25 < z < 0.96 by the first year Supernova Legacy Survey (SNLS) 
and the remaining 119 points are old data. The SN la observation gives the distance modulus of 
a SN at the redshift z. The distance modulus is defined as 



^ithiz; P, M) = 5 logio(dL(^)/Mpc) + 25 = 5 logio [(1 + z) 



dz' 
E(7 



+ 25-51ogioi^o, (7) 



where the luminosity distance (Il^z) = ^^^^p- /g^ "E^' nuisance parameter M = 51ogiQi7o is 
marginalized over by assuming a flat prior P{Hq) = 1 on Hq, P = {c, describes a set 

of parameters characterizing the given model. In order to place constraints on the interacting 
holographic dark energy model, we perform statistics for the model parameter P 



(8) 



Our analysis shows that if we use the SN la data, the constraint is not good, and the la range 
is rather large. 
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(a)l a range of q{z) (b)l a range of w{z) 

Figure 7: (a)The evolution of q(z) within the 1 a range. (fe)The evolution of w{z) within the 1 cr range. The 
1 a range of each combination is between the same colored dash lines. And the solid lines are the best-fit 
curves of each combinations. 



An efficient way to reduce the degeneracies of the cosmo. 



ogical parameters is to use the SN la 



12l ] and the CMB shift parameter 



data in combination with the BAO measurement from SDSS 
The acoustic signatures in the large scale clustering of galaxies yield additional test for cosmology. 
Using a large sample of 46748 luminous, red galaxies covering 3816 square degrees out to a redshift 
of z = 0.47 from the SDSS, Einstein et al 
which is described by the A parameter 



12] have found the model independent BAO measurement 



, , 1 fZBAO rl^l , „ 

A = V^E[zsAo)-"'[^ / 4h^'" = 0.469(-^)-0-35 ± 0.017, (9) 

ZBAO Jo E{z') 0.98 



ZBAO Jo 

where can be taken as 0.95 



22(1 and zbao = 0.35. In our analysis we first investigated the joint 
statistics with the SN la data and the BAO measurement. The result is shown in Figure 1, where 
we show the contours of 68.3% , 95.4% and 99.7% confidence levels. The fitted parameters with 
the la errors are shown in Table 1. 

We also use the CMB shift parameter given by 

"-^-L m-y (^°' 

where zis = 1089. This CMB shift parameter R captures how the /-space positions of the acoustic 
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(a) (b) 
Figure 8: (a)The evolution of entropies with the — 0.08 and c = 1 and the initial conditions ^Ido = 



0.7 and HqSq ~ 10 (6) The evolution of entropies with the best fit parameters of the combination 



Lookbacktime+SN la+BAO+CMB, 6^ = -0.003 and c = 0.83 and the initial conditions floo = 0.71 and 
HiSo = 10-30. 

peaks in the angular power spectrum shift. Its value is expected to be the least model independent 



Hi. 



and can be extracted from the CMB data. The WMAP3 data [22] gives R = 1.70 ± 0.03 
Now we can combine the SN la, WMAP3 and SDSS data to constrain the interacting holographic 
model. Using the statistics, contours from the joint constraints SN la+BAO+CMB are shown 
in Figure 2. Comparing with Figure 1, we see that the errors have been reduced significantly in 
the joint analysis. The la range of the model parameters are listed in Table 1 for comparison. 

It is of interest to include the Hubble parameter data to constrain our model. The Hubble 
parameter depends on the differential age of the universe in terms of the redshift. In contrast to 
standard candle luminosity distances, the Hubble parameter is not integrated over. It persists fine 
structure which is highly degenerated in the luminosity distance [4]. Observed values of H{z) can 
be used to place constraints on the models of the expansion history of the universe by minimizing 
the quantity 

^2^p^ = ^ lHoM-H,,{z,;Pr _ (11) 



This test has been used to constrain several cosmological models 



However this test on its 



own cannot provide tight constraint on the model. It is interesting to combine the H{z) data with 
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(a)SN la + BAO (b)SN la + BAO + CMB 
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ic)H{z) + SN la + BAO (d)Lookback time + SN la + BAO 
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Figure 9: The contours in the holographic dark energy model without interaction. This graph shows 
observational contours in the (c — CId) plane. 
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Table I: The best-fit data of the interacting holographic dark energy model. 
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Table II: The best-fit data of the noninteracting holographic dark energy model. 
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the data above to obtain tighter constraints on the interacting holographic dark energy model. 
The result on the joint analysis H(z)+SN la+BAO is shown in Figure 3 and the la ranges of 
different parameters are listed in Table I. Because the sample of H(z) data is too small at this 
moment, the constraint on the model by including H{z) data is not very tight. We hope that the 
future observations can offer more data of H{z) so that can be reduced. Adding the CMB shift 
parameter data, we have shown the combined analysis H{z)+SN Ia+BAO-|-CMB shift in Figure 
4. Comparing with Figure 3, it is interesting to notice that errors of model parameters have been 
significantly reduced. 



Table III: The best-fit data of ACDM . 
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The constraint based on SN la data and the recently proposed angular-redshift relation of 
compact radio sources are distance based methods to probe cosmological models, now we are 
going to test the model by using the time-dependent observable, the lookback time. The new test 
is expected to provide a complementary test of the model. This method has been employed in 
y, Q, y, [3] • The lookback time -redshift relation is defined by 

where H^^ = 9.78h~^ Gyr, and we use the present value of h = 0.72 given by the HST key project 



23], P stands for the model parameters. To use the lookback time and the age of the universe to 
test a given cosmological model, let's follow j3| to consider an object i whose age ti{z) at redshift z 
is the difference between the age of the universe when it was born at redshift zp and the universe 
age at z, 

/X rr u f°° dz' f°° dz' 

Using the lookback time definition, we have t{zi) = tL{zF) — tiiz)- Thus the lookback time to an 
object at Zi can be expressed as 

tt{zi) = tlizF) - t{Zi) = [tf - U{Z)] - [tf' - tL{zF)\ = t""' - U{Z) - df, (14) 

where df = t°^'^ — tiizp) is the delay factor. 

In order to estimate the parameters of our model, we minimize the function 

Xage\^ ) /-^ ^2 , ^2 _2 ' V-^^ j 

* o o 

where Oi = \ Gyr is the uncertainty in the individual lookback time to the ith galaxy cluster of 
our sample and fi^obs = 1.4 Gyr stands for the uncertainty on the total age of the universe until 
now. The current age of the universe in our analysis is taken as 14.4 Gyr. The second term in 
the expression was introduced to make sure that the cosmological model can estimate the age 
of the universe at present in addition to describing the age of the universe at high redshift. Since 
the delay factor df does not appear explicitly in the theoretical value of tL{zi), we will treat it 
as a nuisance parameter and marginalize it in our calculation. The joint statistical analysis of 
the combined observations including lookback time+SN la+BAO has been done and the result 
is shown in Figure 5. Comparing to the analysis of SN la+BAO shown in Figure 1, we noticed 
that the parameter space now is enlarged. This fact is expected and understood in term of the 
conservative uncertainty assumed (cii = IGyr) for the individual lookback time. In Figure 6, we 
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have shown the combined analysis including lookback time+SN la+BAO+CMB shift. It is easy 
to see that adding the CMB shift data, the model parameters have been constrained much tighter. 

To illustrate the cosmological consequences led by the observational constraints, we show the 
evolution cases of the equation of state parameter w{z) and the deceleration parameter q(z) accord- 
ing to the best-fit values of our model parameters in Figure 7. It is easy to see that our model can 
have the feature of w crossing —1. Our present equation of state and the deceleration parameter 
are consistent with CMB data [6, 221. 



It is interesting to notice that our best fit value of 5^, the coupling between dark energy and 
dark matter, is negative. In the holographic interacting dark energy model by employing the 
apparent horizon as the IR cutoff [2^, it was argued that an equation of state of dark energy 
tt; < is necessarily accompanied by the decay of the dark energy component into pressureless 
matter (6^ > 0). However, in our model, negative 6^ can accommodate reasonable equation of 
state of dark energy which is clearly shown in the Figure 7. Another worry of the negative b'^ 
which implies a transfer of energy from the matter to the dark energy is that it might violate the 



second law of thermod yna mics 
employ the formula in 



2ll |. In order to check the second law of thermodynamics, we can 



171 ]. Using the apparent horizon as a thermal boundary and evaluating the 
entropy inside the apparent horizon from the Gibbs law, we have shown the evolution of entropies 
in Figure 8. It is easy to see that for the best fit negative 6^, entropy of matter and fluids inside 
the apparent horizon plus the entropy of the apparent horizon do not decrease with time. The 
generalized second law of thermodynamics is still respected. 

For the sake of comparison, we have shown the same contours of the holographic dark energy 
model without interaction in Figure 9. And the best-fit results are shown in Table II. Combined 
SN Ia-|-CMB-|-BAO constraints on the holographic dark energy model without interaction have 
been studied in ^2^]. Here we have added the Hubble diagram data which is not an integrated 
over effect and the time-dependent observable analyses. Comparing with Figure 2, we find that 
2a, 3a confidence ranges in the c — Qdo plane are much smaller for the holographic dark energy 
model without interaction. In the analysis we find that the difference between the model with 
and without interaction is bigger when the model parameter c is bigger. We have also listed the 
best-fit results for a flat ACDM model in Table III. At the flrst glance, from the Xmim gives 
us a sense that the interaction between dark energy and dark matter gives a better description of 
the combined observations although this interaction is extremely small. Considering the additional 
degree of freedom, it is still early to say that our interacting holographic dark energy model is more 
favored than the ACDM model. However one advantage of the interacting holographic model is 
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that, unlike the ACDM model, it can alleviate the coincidence problem 

In summary, in this work we have performed a parameter estimation of the interacting holo- 
graphic dark energy model which could explain the observed acceleration of our universe. We have 
analyzed data coming from the most recent SN la samples, CMB shift, LSS observation, H{z) 
and lookback time measurements. Comparing with the single observational test, we learnt that 
the joint analysis of different observations based on different observables is powerful to overcome 
the statistical uncertainties. We have got useful consistent check of the interacting holographic 
dark energy model and tighter constraints on the model parameters. The joint analysis indicates 
that this is a viable model. In the Icr range, it can explain the transition of the equation of state 
from w > —1 to w < —1. It is worth noting that although the current H{z) and lookback time 
data do not provide very restrictive constraints, richer samples of H(z) data and more precise 
age measurements of high-z objects will provide a complementary check of the cosmic acceleration 
model. The joint statistical analysis is necessary to test the model. 

We observed that the best fit coupling between dark energy and dark matter is negative, which 
indicates that there is a possible energy transfer from the matter to the dark energy. Although the 
generalized second law of thermodynamics is shown not threatened by the best fit negative 5^, the 
holographic principle might still violate in the future if there is a continuing energy transfer from 
the matter to the dark energy, since in the late stage it is possible to see that S < Sa in Figure 8. In 
view of the unknown nature of dark energy and dark matter, we can't say for certain the direction 
of the energy transfer between dark energy and dark matter. However from the observation data 
and the holographic principle requirement, the nature thing we can learn is that in the past there is 
an energy flow from the matter to the dark energy from the observational data, while in the future 
there requires an energy transfer from the dark energy to the matter to satisfy the holographic 
principle. A nature description of the coupling between dark energy and dark matter is called 
for, which is important to influence the structure formation and the description of the universe 
evolution. 
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